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Abstract: We investigate semiclassical properties of space-time geometry of the low en¬ 
ergy limit of reduced four dimensional supersymmetric Yang-Mills integrals using Monte 
Carlo simulations. The limit is obtained by a one-loop approximation of the original 
Yang-Mills integrals leading to an effective model of branched polymers. We numerically 
determine the behaviour of the gyration radius, the two-point correlation function and the 
Polyakov-line operator in the effective model and discuss the results in the context of the 
large-distance behaviour of the original matrix model. 
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1. Introduction 

The ten dimensional matrix model of reduced supersymmetric Yang-Mills integrals [1] is 
believed to be a good candidate for a non-perturbative definition of string theory for the 
following reasons [2, 3, 4], In the large N limit the action of the model takes the form 
of the action of IIB strings. The model is supposed to contain in the matrix structure all 
topological excitations of the string world-sheet. The one-loop approximation of the model 
suggests that the space-time geometry dynamically undergoes a spontaneous dimensional 
reduction from ten to four dimensions - a spontaneous breaking of the ten dimensional 
Lorentz symmetry. It was actually the first ever proposed candidate in string theory for a 
dynamical mechanism of dimensional reduction, so the model has attracted the attention 
of many researchers. In the context of string theory it is called IIB matrix model or IKKT 
model after the authors [1]. 

The conjecture about the spontaneous symmetry breaking has been derived in the one- 
loop approximation of the model [3, 4], This approximation leads in the low energy limit 
to an effective model of graphs which are dressed with vector fields representing positions 
of string world-sheet points in space-time. The graphs are closely related to branched 
polymers which have fractal dimension four [6, 7, 8]. The essence of the conjecture is that 
starting with a ten dimensional model one ends up with a four dimensional one [3, 4]. 

It has become crucial to check whether the one-loop approximation is a good approxi¬ 
mation in the limit of large distances and whether the scenario holds beyond the one-loop 
approximation [9, 10, 11, 12]. In this paper we study the first issue. 

The model is not solvable analytically and it is also very complicated to treat nu¬ 
merically because of a sign problem which appears when one integrates out the fermionic 
degrees of freedom [13]. Many approximation techniques have been developed to investi¬ 
gate the behaviour of the model [3, 4, 5, 14, 15, 17, 18]. Because of the complexity of the 
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ten dimensional model its four and six dimensional counterparts have been studied to gain 
deeper insight into the geometrical nature of the problem and the limitations of effective 
low energy models [16, 17, 18, 19]. The four dimensional case is much easier since it is 
free of the sign problem. In four dimensions one has discovered some universal properties 
of the model which also hold in higher dimensions. For instance one has found a relation 
between the existence of elongated spiky world-sheet configurations and the singularities 
of the partition function. A simple argument about entropy of spiky geometries derived in 
four dimensions could be extended to higher dimensions. It has allowed for deducing the 
proper singularity type of the partition function, unravelling the geometrical structure of 
configurations responsible for the singularities [18]. 

In this paper we also concentrate on the four dimensional case to analyse the geomet¬ 
rical behaviour of the low energy effective model given by a model of branched polymers 
obtained from the supersymmetric Yang-Mills matrix model in the one-loop approxima¬ 
tion. We numerically study the gyration radius, the two-point correlation function and the 
Polyakov-line operator whose effective form was derived analytically in [20] . 

Four dimensional supersymmetric Yang-Mills integrals are interesting per se and in 
connection with supersymmetric QCD. Actually the relation of reduced integrals to the 
1/lV-expansion was first discovered in QCD [21] where it was also applied to study the zero 
momentum content of the theory [22]. 

2. The model 

The partition function of the model is given by a reduced supersymmetric Yang-Mills 
integral [1, 2, 3] 

Z = j ZL4DTDT e ~ s[A ^^ ] , (2.1) 

with the action 

S , [A,'F,'P] = -j^Tr [A^A"] 2 - -^TV r“ b [A^,'I' 6 ] , (2.2) 

where A (// = 1,... D) are traceless N x N Hermitian matrices and \k a are N x N 
traceless matrices of Grassmannian variables, transforming as Majorana-Weyl spinors for 
D = 10, or as Weyl spinors for D = 4. We shall discuss here the D = 4 case. 

One is interested in the large N behaviour of the model. It can be studied by approxi¬ 
mate methods. In particular one can use perturbation theory to estimate the contribution 
of quantum fluctuations around classical solutions [3, 4], The main idea consists in splitting 
the fields yp 1 , 'P® and \P a into classical part and quantum fluctuations as 

AC = xfSij + aC % = tfdij + r't + r" , (2.3) 

and integrating out the quantum fluctuations a and ip. In addition, it is also convenient 
to integrate out the classical fermionic fields £ to get rid of the Grassmannian content in 
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the effective model. Having done this one is led to a model of branched polymers with the 
partition function 

Z bp= E f Hd 4 Xi e~V’ x] , (2.4) 

T£T n J i 

which describes the large-distance behaviour of the semiclassical fields xf (we shall write 
in short Xi) representing the positions of string world-sheet points in the target space. If 
the positions are widely separated, that is |x) — Xj\ 3> y/g, the effective action is given by 
[3,4] 

5 bp [T,x]=6 E ln\xi~Xj\ • (2-5) 

A detailed analysis of the matrix model shows that a strong repulsion occurs whenever any 
two vectors Xi and Xj come close to each other [3, 4], One believes that details describing 
the short-distance repulsion do not affect the universal large-distance properties of the 
system. One can therefore choose the way in which one models the repulsion as long as it 
is short-ranged. One can for instance use a hard core repulsion, which is relatively easy to 
implement. In this case the action (2.5) reads 

{ 6 ln|xj — Xj\ if \x a — Xb\ > c Va,b€T 

(H)eT ( 2 . 6 ) 

+oo otherwise 

The size of the core c is a free parameter in the model. If any two vertices of the branched 
polymer come too close to each other the action becomes infinite and the corresponding 
configuration is entirely suppressed in the partition function. Each vertex is surrounded 
by a ball-shaped zone which contains no other vertices of the branched polymer. From 
the numerical point of view the branched polymer model is much easier to simulate than 
the matrix model. The complexity of the algorithm, defined as the number of operations 
to update all degrees of freedom grows as IV 8 for the SU(N ) supersymmetric Yang-Mills 
model (2.1) and as N 2 for the branched polymer model (2.4). In the first case the number 
of degrees of freedom which one has to update in one sweep is proportional to the number 
of elements of the A^ matrices, and grows as N 2 . In an update of an element of the matrix 
one has to compute a determinant of the Dirac operator which is a K x K matrix of size 
K oc N 2 . Computation of a determinant requires oc K 3 operations, yielding altogether 
an at least IV 8 complexity of a sweep in the matrix model. In the second case in a sweep 
one has to update N vertices. In each update one has to make N — 1 global checks of the 
hard core condition with remaining vertices. This yields an N 2 complexity of a sweep in 
the branched polymer model. The presence of a hard core leads to an additional slowing 
down of the algorithm: Branched polymers have generically fractal dimension four. Thus 
if they are embedded in four dimensions they densely fill up the space. Because of the 
hard core constraints the configuration space looks like a cheese with empty holes and the 
algorithm has to maneuver between them while constructing new configurations. Since the 
algorithm does it by trial and error it takes a long time to produce completely uncorrelated 
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configurations. Even if one takes the autocorrelations into account the complexity of the 
branched polymer algorithm is many orders of magnitude smaller than of the corresponding 
matrix model algorithm for the same N, already when N is of order ten. 


3. Physical quantities 


Physical quantities are defined as operators 0[^4 M , \k] which depend on the matrices Ap and 
T in the original model (2.1). In the one-loop approximation, after integrating out quantum 
fluctuations a M , ^ and classical fermionic fields £, the effective operators O\A ,X , vp] O bp [x\ 
become functions of a branched polymer graph T and fields x t dressing its vertices, which 
we shall denote in short by x. For each operator one has to derive its branched polymer 
counterpart. The Polyakov-line operator, which is the fundamental operator in the theory 

Pk = -^Trexp [ik^A**] (3.1) 

is mapped into the following operator in the branched polymer picture [20]: 


Pbp,k [%\ 


1 

N 




\ ^ \ ^ _ 

a ik p( X “ ~ X P b )\x a - X b \ 2 


(3.2) 


e ik^x^ __ e ik„x% 


( c d^T ik P^ ~ X 0K ~ x b \ 2 


21 xg - x b \ 2 k 2 - ((x% - xp kg-Y 
3|x a - x b \ 2 


The operator contains three terms. The first one has a classical origin. It is given by a 
sum over branched polymer vertices. The second contains a sum over all pairs of vertices 
independently of whether they are neighbours on the branched polymer or not, while the 
third one - a sum over pairs of neighbouring vertices on the branched polymer graph. These 
two terms come from quantum fluctuations in the one-loop approximation. 

The most fundamental physical observable describing the distribution of the classical 
fields Xi in the target space is the two-point correlation function: 


cf (W X 2 ) = / 2 ]T d(xi - X{) ]T S(xj -X 2 )\ . (3.3) 

\ ^ i<3 j / N 

The average is taken over the ensemble {T, Xj} of branched polymers with N vertices with 
the partition function (2.4). The statistical weight of the configuration is given by a product 
of link weights 

e~ S = n Wdxi-Xjl) (3.4) 

(H)eT 

on which additionally the hard core constraints are imposed. Since the statistical weight 
and the hard core constraints are invariant with respect to a shift of all coordinates by the 
same vector X{ —> Xi+5, the two-point function G y N depends on the difference X = X 2 —Xi, 
or actually - due to the isotropy of the statistical weight - only on the distance X = |X|. 
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Using the two-point function one can determine the typical linear extent of the system. 
Usually one does it by calculating the second moment of the two-point function: 

X)« = jT dXX 2 C$(X) = ^ JV(JV 2 _ D Y, Pi - %l 2 y ■ (3.5) 

The square root of this expression gives a quantity called radius of gyration which is a 
standard measure of the linear dimension. For branched polymers one can change integra¬ 
tion variables in the partition function (2.4) from x's to r’s: 'f)y = x t — Xj. Now one can see 
that the integration over r^-’s in (3.5) leads to a divergence when it is done for link weights 
with power-law tails W(r) ~ r -6 : the integration of r 2 over the four dimensional volume 
d 4 r = fi r 3 dr, where Q is the angular part of the integration measure, gives a logarithmi¬ 
cally divergent quantity: (r 2 ) = U f drr 3 r 2 W(r ) ~ f dr/r. The gyration radius (3.5) is ill 
defined. In this case the linear extent can be defined by the first moment of the two-point 
correlation function [8, 17]: 

Rn = (X) N = dXXG%\x) = ^ jy(jV 2 _ 1} J2 ■ ( 3 - 6 ) 

We will use this quantity in this paper to measure the linear extent of the system. 

4. Geometry of branched polymers 

We generate our configurations using a Monte Carlo algorithm. It consists of two update 

types which are applied alternately: a graph structure update and a vertex positions update 

[8]. When the branched polymer topology is updated vertex positions are kept constant, 

while when the vertex positions are updated the graph’s topology is kept constant. We 

2 

tested the algorithm for Gaussian branched polymers with link weights W(r) = e~ r , 
r = \xi — Xj\ and without a core. In this case one can determine the form of the two- 
point function analytically [8]. For large N one expects the two-point function G‘ N (X) to 
effectively be a function of the scaling variable X/N 1//4 . In other words one expects that the 
fractal dimension is Df = 4. One can see in figure 1 that indeed the data for the two-point 
functions for different N collapse to one curve if one rescales X —> X/N 1 / 4 . We now turn to 
the branched polymers (2.6) derived from the matrix model (2.1). The two-point function 
for different system sizes is plotted against the scaling variable in figure 2. One can see 
small deviations from the scaling. The reason for this behaviour is related to the presence 
of the power-law tail in the link length distribution W (r) which leads to the divergence of 
the second moment of the distribution (r 2 ) = f d 4 r r 2 W(r) ~ f dr/r as already mentioned 
above. This means that fluctuations of the link lengths are infinite and that from time to 
time long links appear in the system. These links are reminiscent of the spiky configurations 
observed in the surface model [18]. Let us now try to understand algorithmical problems 
caused by the presence of such a fat-tailed distribution. The algorithm which we use is 
an example of dynamical Monte Carlo. It generates a Markov chain of configurations 
randomly walking in the configuration space. Any two consecutive configurations in the 
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Figure 1: The two-point correlation function G^\x) for the Gaussian branched polymer model 
plotted against X/N 1//4 . 


Figure 2: 



power-tail weights and a hard core c = 1 plotted against X/N 1 / 4 . 


chain differ from each other only a little. In a single step the algorithm changes a link 
vector by a small value r —> r + 5 and accepts this change with the Metropolis probability. 
The sequence of changes can be viewed as a sort of random walk in the potential In W(r). 
This potential is very flat and therefore if the algorithm once produces a long link it takes 
a long time to make it short again. Such a random walk algorithm introduces therefore 
large autocorrelation times. Moreover, a long excursion towards the tail of the distribution 
produces many long links, so after each long excursion a surplus of long links and a deficit 
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of short ones occurs in the recorded history. This is an effect which makes the rescaled 
histograms in figure 2 not to coincide. The narrow tails of the histogram go far beyond 
the range in the figure. The tail behaviour of the two-point function is presented in the 
logarithmic scale in figure 3. It is compared with the tail behaviour of the original matrix 
model (2.1). The two-point correlation function inherits the power-law tail drr ~ 3 from 
the link length distribution ~ drr^W(r). As we can see in figure 3 in both cases the data 
behave as 



power-tail weights regularised by introducing a hard core c = 1 and the IIB matrix model for 
different values of N. The solid black line illustrates the expected power-law behaviour: ~ X~ 3 . 

The presence of the tails ~ r -3 makes the determination of the linear extent of the 
system R jy (3.6) difficult. The computation of Rjy amounts to the computation of the 
first moment of the two-point function. The autocorrelation time for the algorithm is large 
and grows with N. On the other hand one expects using the central limit theorem that 
the tail r ~ 3 belongs to the Gaussian universality class, what means that the probability of 
entering the tail part of the distribution decreases with the number of degrees of freedom, 
in our case with N. One therefore expects that in the limit N —> oo the quantity Rn 
should depend on the bulk of the distribution and not on the tail which becomes marginal 
in this limit, and thus that the broadening of the statistical error coming from the large 
autocorrelation time should be finite for Rjy for large system sizes. One can use standard 
methods to estimate the error bars of Rjy. If we do so we obtain the data presented in 
figure 4 illustrating the dependence of the linear extent Rn on the system size. The data 
presented in figure 4 can be fitted with the scaling formula fi(N) = ai • IV 1 / 4 . If we apply 
this formula to find the best fit for the range N > 48 we obtain ai = 3.44(5) with the fit 
quality: Xdof = 2-9- If we include finite size corrections / 2 (IV) = <22 • IV 1 / 4 (1 + 62 IV' 1 ), 
then the best fit is for 02 = 3.56(4) and 62 = —6.8(9) and has the quality Xdof = 11- 
The fit quality improves in comparison with the previous one as can also be seen with the 
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Figure 4: Comparison of the average gyration radius Rn for the branched polymer model with 
a hard core c = 1 plotted against the number of vertices N. The continuous lines represent the 
functions: fi(N) = 3.44 • N 1//4 and ^(IV) = 3.56 • iV 1/,4 (l — 6.81V -1 ). 

bare eye. The finite size corrections which control deviations from the pure scaling can 
be attributed to the power-law tails which now and then cause the appearance of a long 
link on the branched polymer and to the hard core repulsion which also effectively makes 
the polymer look more elongated. The repulsion disfavors crumpled trees. In principle one 
could expect the fractal dimension of branched polymers to lower. The magnitude of the 
finite size corrections in / 2 (IV) is for bo = — 6.8 of order a few percent for N is of order a few 
hundred. One can thus believe that the results are consistent with the fractal dimension 
Df = 4. To thoroughly check this hypothesis one would though have to study much larger 
systems. 

5. Polyakov-line operator in the branched polymer model 

In this section we will discuss measurements of the Polyakov-line operator. As it stands 
in equation (3.2) it has real and imaginary parts. The imaginary part however gives zero 
on average due to the reflection symmetry: the partition function of the branched polymer 
model is invariant with respect to simultanous reflections Xi —> —Xi of all vertex positions 
and therefore the average value of the Polyakov-line operator over the ensemble of branched 
polymers fulfills the condition: 




(5.1) 


Thus it is sufficient to measure only the real part of the operator. We shall denote it by 
Pf.. The momentum k plays the role of an external parameter. Since the branched polymer 
system is isotropic one can choose k to lie on one of the axes of the coordinate system in 









which components of the fields Xi are expressed. In practice, while doing numerical simu¬ 
lations one can calculate projections on the four independent components // = 1... 4 and 
average over them to improve statistics. This procedure eventually leads to the following 
operator 

Pk M = M + Pl 2) [x] + if > [x] (5.2) 


where 

p fc (1) w = ’ 

/i=l a 

p (2) r 1 _ -k -A sin (kx%) - sin(kx£) 

k x 4IV _ x ^)\x a - Xb \ 2 

/x = 1 a<b x 0/1 1 

( 2 ) \ , _ sin(fcxq) - sin(fc.T^) 2|x a - x b | 2 - (a# - x^)~ 

X { x a -X%)\x a -X b \ 2 3|x a - X b \ 2 

In the last equations k stands for the length of the vector k, which is independent of the 
space-time direction fi. In particular kx M is a product of the length of the vector k and the 
/j-th component of x. 

The length scale in the branched polymer model is set by the size of the core c. It is the 
only scale parameter in the model. The model is invariant under a simultaneous rescaling 
c —> Ac and x,{ —» Ax). One can use this invariance to rescale the core size to c = 1. The 
average of an operator 0[x\ over the ensemble of branched polymers with a hard core of 
size c can be expressed as the average of this operator for the rescaled argument 0\cx\ over 
the system with a core c = 1: 


{0[x\)c = {0[cx\)c= 1 . (5.3) 

In particular for an operator which is a homogeneous function of order A: O [Ax'] = A a O[x] 
we have 

(0[x]) c = c A (0[x]) c= i. (5.4) 

The momentum k of the Polyakov-line operator introduces an additional scale which has 
to be measured relatively to c _1 . If one rescales x Ax one should simultaneously rescale 
k —» k/A to keep the combination kx constant. We can now relate the value of the 
Polyakov-line in the system with a hard core c to its value in the system with c = 1: 



The amplitudes of the second if [x\ and the third term pjf' 1 \x\ scale as c 4 with the core 
size while the amplitude of the first term pf [x] stays constant. Using this observation 
we can reconstruct the behaviour of the average of the Polyakov-line operator (Pk)c in the 
system with an arbitrary c from measurements of (pf), (pf) and ( P'^) for c = 1. We 
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Figure 5: The three terms {Pj. 1 ' > )c=i, (P^ 2 ^) c =i and (P^) c = 1 contributing to the average Polyakov- 
line operator in the one-loop approximation (5.2) in the system with hard core size c = 1. 


will follow this strategy. In figure 5 we present the dependence of the three contributions 
for c = 1 on the momentum k, obtained by numerical Monte Carlo measurements. Since 
the second and the third terms are negative, the Polyakov-line operator being a linear 
combination of the three terms (5.2) may assume negative values for some range of k 
depending on the value of c (5.5). Such a behaviour is also observed in the matrix model 
for small N as we can see in figure 6. In the matrix model the effect however disappears 
when N increases and is practically absent already for N of order ten. The curve for N = 2 
in figure 6 was obtained analytically as a Fourier transform of the eigenvalue distribution 
whose form is known for SU( 2) [23]. 

Now the issue is to test up to what degree the model of branched polymers can re¬ 
produce the full matrix model. The model of branched polymers was constructed as an 
effective model for the large-distance behaviour of the original model. So we are now look¬ 
ing for the optimal value of the parameter c, which is a free parameter of the effective 
model, to reproduce the behaviour of the matrix model on large distances which corre¬ 
spond to small values of k in the Polyakov-line operator. For the comparison we will use 
the matrix model results presented in [17]. They correspond to N = 16, 24, 32 and 48. 
The coupling constant g of the matrix model is chosen to be g = y/ 48/IV and the results 
are expressed in physical units k phys = k/y/g [17]. After rescaling to physical units the data 
of the matrix model for different N collapse to one curve. This curve goes as 1 — ak 2 + ... 
for small k. We found that for c = 0.75 also the branched polymer results follow this curve 
in the region of small momenta. The results are presented in figure 7. Indeed for small 
k the data points for different N follow the same line 1 — ak 2 . The slope parameter a is 
thus universal. For larger k 2 the results depart from the universal line in a non-universal 
manner which depends on N. There are two reasons for that. First of all the effective 
model by construction is only supposed to describe the large-distance (small k) behaviour. 


- 10 - 









Figure 6: Comparing the Polyakov one-point correlation functions (Pk) of the IIB matrix model 
for N = 2, 3, 4, 5, 6. 



k 2 /g 


Figure 7: Comparing the 1 — ak 2 + ... behaviour of the Polyakov-line operator (Pk) of a branched 
polymer with c = 0.75 for different sizes of N and the IIB matrix model results. 


For larger k the effective model depends on the details of the regularization and is not any 
more universal. Secondly the Polyakov-line operator (3.2) contains only one-loop correc¬ 
tions. Higher loop corrections, which have been neglected in (3.2) would introduce terms 
of order 0(fc 4 ). Let us stress again, what is most important in figure 7 is that the leading 
k 2 behaviour is universal and identical with that of the original matrix model. 
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6. Summary 

We have shown that the branched polymer model very well describes all essential features 
of the reduced supersymmetric Yang-Mills integrals in four dimensions at large distances: 
the power-law tail behaviour of the two-point function, the four dimensional scaling of 
the gyration radius and the universal small momentum behaviour of the Polyakov-line 
operator. The branched polymer model can thus be useful in studies of the low momentum 
content of supersymmetric QCD [22], 

What remains to be checked is whether the one-loop approximation works equally well 
for the ten dimensional matrix model where branched polymers are replaced by somewhat 
more complex graphs in the effective model. These graphs have branched polymer back¬ 
bones and one believes that they preserve the four dimensional fractal structure. Additional 
links which appear on the graphs introduce some rigidity constraints which therefore may 
make the fractal dimension to become the canonical dimension of the underlying geome¬ 
try. This would explain the existence of a dynamical spontaneous symmetry breaking of 
the original ten dimensional Lorentz symmetry. There are some indications that such a 
scenario indeed takes place [9, 10, 11, 12, 13, 24, 25]. 
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